Report for Model 10

— for both jointly and separately

2024-08-19

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.03
## LDevsum       1.01       1.05
## dh0           1.38       2.18
## dh1           1.72       3.73
## dh2           1.47       2.65
## dl0           1.09       1.34
## dl1           1.19       1.68
## dl2           1.03       1.03
## dl3           1.19       1.66
## 
## Multivariate psrf
## 
## 1.38

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1227.65640 22087.6761
DIC3 1160.16296 22151.3761
PWAIC 44.21523 242.1802
WAIC 1187.89048 22173.1181

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.08       1.33
## dh0           1.41       2.23
## dh1           2.51       5.00
## dh2           1.57       3.22
## dl0           3.28       6.91
## dl1           1.07       1.08
## dl2           3.35       7.44
## dl3           1.76       3.22
## 
## Multivariate psrf
## 
## 3.05

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1381.02334 22212.1649
DIC3 1246.20600 22234.6104
PWAIC 90.86349 297.1184
WAIC 1315.67881 22266.9946